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Abstract 

The subscale (11- and 24-inch) hybrid 
motore at the Marshall Space Flight Center (MSFC) 
have been used as versatile and cost effective testbeds 
.for developing new technology, • Comparisons 
between motor configurations, ignition systems, feed 
systems, fuel formulations, and nozzle materials have 
been carried out without detailed consideration as to 
how "good” the motor performance data were. For 
the 250,000 Ib r -thrust motor developed by the Hybrid 
Propulsion Demonstration Program consortium, this 
shortcoming is particularly risky because motor 
performance will likely be used as part of a set of 
downselept criteria to choose between competing 
ignition and feed systems under development. 

This analysis directly addresses that , 
shortcoming by applying uncertainty analysis 
techniques to the experimental determination of the 
characteristic velocity, theoretical characteristic 
velocity, and characteristic velocity efficiency for a 
24-inch motor firing. With the adoption of fuel-lined 
headends, flow restrictors, and aft mixing chambers, 
state of the art 24-inch hybrid motors have become 
very efficient However, impossibly high combustion 
efficiencies (some computed as high as 108%) have 
been measured in some tests with 11-inch motors. 
This analysis has given new insight into explaining 
how these efficiencies were measured to be so high, 
and into which experimental measurements contribute 
the most to the overall uncertainty. 

English 

a: Initial throat diameter for erosion model 


A*; Nozzle throat area (in 5 ) 

b: Fitting constant for erosion model 

B,: Bias uncertainty 

Bi; Biot number 

C*: Characteristic velocity (fi/sec) 

C d : Discharge coefficient 

D: Diameter (in) 

DRE: Data reduction equation 

OOX: Gaseous oxygen 

GUI: Graphical user interface 

h: Perturbation magnitude 

ty Hot gas-side convective heat transfer 

coefficient 

HFDF: Hybrid Propulsion Demonstration Program 
I ff : Sea-level specific impulse (lbf-sec/lbm) 

I Vie : Vacuum specific impulse (lbf-sec/lbm) 

JIRAD; Joint Independent Research and 
Development 

k Thermal conductivity (^TU/in-sec-Tl ) 
LOX: Liquid oxygen 

LSSRCS: Large Subscale Solid Rocket Motor 
Combustion Simulator 
M: Mass (lbm) 

m : Mats flow rate (Ibm/iec) 

NUj,: . Nuaselt number based upon diameter 
ODE: One-dimensional equilibrium 

O/F: Oxidizcr-to-fud ratio 

P: Pressure (ftf/in 1 ) 

Pr: Frandtl number 

P,: Precision uncertainty 

r Experimental result 

Rep: Reynolds number based upon diameter . 

SHUE: Software for Hybrid Uncertainty Estimation 

SLOC: Source lines of code 

t*,: Motorbumtime 

TEA/TEB: Pyrophoric mixture of 

triethylaluminum and triethylboride 


* Propulsion Engineer, Advanced Space Systems, Member AIAA 

* Eminent Scholar in Propulsion and Professor, Associate Fellow AIAA 
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UPC; Uncertainty percent contribution 

U,: Experimental uncertainty 

X<: Primitive measured variable 

Greek symbols 

a: Coefficient of thermal expansion (l/ 6 ^) 

AT: Temperature difference (°R) 

<j>: Sum of 0, 0 3 mass fractions in exhaust 

species 

Y Ratio of specific heats 

T)^: Combustion efficiency 

\l: Viscosity (Ibm/iii-aec) 

p: Density (lbm/ft ! ) 

Subscripts and Superscripts 

*: Nozzle Throat 

Vap: Vapor 

0: Chamber condition 

laflaMiaa 

Lockheed Martin Astronautics’ (as the 
former General Dynamics, then Martin Marietta, 
before merging to form the ament Corporation) 
involvement with hybrid propulsion development at 
the Marshall Space Flight Center dates back to 1990, 
with the formulation and execution of the JIRAD 
program. 

The successful JIRAD program effort was 
followed in 1992 by the LSSRCS program, a concept 
that used subscale hybrid motors to simulate the 
combustion process and exhaust species of a solid 
rocket motor. Numerous test series were completed 
that validated this cost effective evaluation process. 
LMA also supported the HPTLVB program, which 
started in 1993. 

Several hybrid propulsion development 
programs were integrated in 1994 to form a single, 
focused effort to develop and test a series of 250, 000- 
pound thrust hybrid motors (250K configuration). 

Common to all these programs are three 
subscale motors at the NASA Marshall Space Flight 
Center MSFC). Two eleven-inch diameter motor 
configurations (11” case O.D.) are used, one for QOX 
and the other for LOX. . A twenty four-inch 
configuration (24” case O.D.) that exclusively uses 
LOX is also available. The subscale motor 
configurations have matured throughout the different 
programs, but all currently use a headend 
vaporization chamber (lined with either fhel or silica, 
phenolic) and an aft end mixing chamber (similarly 


lined with fuel or silica phenolic, and in the 1 1-inch 
GOX motor, a silica phenolic baffle is also inserted) 
to enhance combustion efficiency and stability. This 
pro- and post-chamber configuration has also carried 
over into the 250K motor. 

.Another common theme among all the 
programs, including HPDP, is that no detailed 
analysis has been done to address measurement 
uncertainties. The fundamental question that 
uncertainty analysis attempts to answer is, “How 
good are the data and the results calculated from 
them?” The actual process is straightforward, but the 
final product requires care to be sure all the sources 
of uncertainty have been accounted for. 

In this effort, a FORTRAN 90 computer 
code was written to calculate uncertainties for 28 
different experimental results for the four hybrid 
motors. The code uses standard analysis techniques 
to estimate partial derivatives and to propagate those 
derivatives, along with their associated measurement 
bias and precision limits, into an overall experimental 
uncertainty. This overall process is demonstrated 
using hot-fire data from a recent HPDP 24-inch motor 
test, 

Additionally, some recent 1 1-inch motor test 
results show impossibly high combustion efficiencies. 
It is asserted in this paper that these impossible 
measurements result from conceptual (modeling) 
biases in the calculation procedures. 

Uncertainty Analysis 

Only a brief overview of the methodology to 
obtain uncertainty estimates and how they propagate 
through a given DRE is given here. The reader is 
referred to Coleman and Steele u for a detailed 
discussion of uncertainty analysis techniques. 

The word accuracy is generally used to 
indicate the relative closeness of agreement between 
an experimentally determined value of a quantity and 
its true value. Error is the difference between the 
experimentally determined value and the truth; 
therefore, as error decreases, accuracy is said to 
increase. Only in rare instances is the true value of a 
quantity known. Thus, it is necessary to estimate 
error, and that estimate is called an uncertainty, U. 
Uncertainty estimates are made at some confidence 
level - a 95% confidence estimate, for example, 
means that the true value of the quantity is expected 
to be within the ±U interval about the experimentally 
determined value 95 times out of 100. 

Total error can be considered to be 
composed of two components; a precision (random) 
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component, 6, and a bias (systematic) component, (}. 
An error is classified as precision if it contributes to 
the scatter of the data; otherwise, it is a bias error; As 
an estimator of (3, a systematic uncertainty or bias 
limit, B, is defined. A 95% confidence estimate is 
interpreted as the experimenter being 95% confident 
that the true value of the bias error, if known, would 
fall within ±B, A useful approach to estimating the 
magnitude of a bias error is to assume that the bias 
error for a given case is a single realization drawn 
from some statistical parent distribution of possible 
bias errors. As an estimator of the magnitude of the 
precision errors, a precision uncertainty or precision 
limit, P, for a single reading is defined, A 95% 
confidence estimate of P is interpreted to mean that 
the ±P interval about the single reading of X, should 
cover the (biased) parent population mean, Ji, 95 
times out of 100. 

In nearly all experiments, the measured 
values of different variables are combined using a 
DRE to form some desired result A general 
representation of a DRE is 

r = r(X l ,X l X,) (1) 

where r is the experimental result determined from J 
measured variables X,. Each of the measured 
variables contains bias errors and precision errors. 
These errors in the measured values then propagate 
through the DRE, thereby generating the bias and 
precision errors in the experimental result, r. 

If the “large sample assumption" is made 1 , 
then the 95% confidence expression for U t becomes 

U) - B) + P? (2) 

where the systematic uncertainty (bias limit) of the 
result is defined as 

(3) 

M M *-/♦! 


and where B ik , the 95% confidence estimate of the 
covariance appropriate for the bias errors in X, and 
is determined from 


(?) 

«■! 

where variables X, and X k share L identical error 
sources. These terms account for correlation between 
bias errors in different measurements. 

For the analysis summarized in this paper, it 
was assumed that bias uncertainties dominated the 
overall measurement uncertainty, making precision 
uncertainties negligible. The only measurements 
where precision uncertainties are important are for 
nozzle throat diameter measurements and the fuel 
weights. As will be discussed, significant non- 
symmetric erosion typically occurs at the throat 
region. These random (from a pre-test prediction 
standpoint) variances in the throat diameter are 
treated mathematically as precision uncertainty 
sources, but in reality, are a unique class of bias 
uncertainty sources termed fossilized biases. This 
mathematical treatment can be found in any statistics 
text, but the 95% confidence large sample precision 
limit is estimated as 




P- =2 

Xi 


Jn 


( 8 ) 


Inspection of equation (3) shows a 
significant effect that correlated bias uncertainties can 
have. Because of the covariance expression of the 
latter half of the right hand side, and most especially 
because the individual partial derivatives are used (as 
opposed to their squares), the correlated bias terms 
can either increase the overall uncertainty (if the 
partials are of the same sign) or decrease the overall 
uncertainty (If the partials are of opposite sign). This 
fact piays a distinct role in this analysis, 


with the bias limit estimate for each X, variable as the Pprffr fjfl flnce Calculations 

root sum square combination of its elemental 

systematic uncertainties Measured C* 



C* is one measure of the energy content in a 
(4) chemical propellant combination. It is defined as 

PA' 

~ (9) 
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However, instantaneous fuel flow rates in 
hybrid rocket motors are currently difficult, if not 
impossible, to measure with any reasonable certainty, 
Therefore, fuel usage is measured by weighing the 
fuel grains (for the 1 1-inch) or the loaded motor case 
(for the 24-inch and 25QK motors) both pre- and post- 
test, and calculating an average fuel flow rate. The 
result is that C* measurements must be based upon 
expended values for fuel, ignitor and inert propellant 
and LOX weights, average nozzle throat area, and a 
total pressure Integral. In effect, all the variables are 
integrated over the motor bum time, resulting in a 
modified from of (9): 


A’ 

’-.,1 \ Pd ' t 


trM 

1 

j M LOX dt + 
J, i 

fro** dt + j jjitf, dt+ j 

>. ». i. 

\mm dt 

0 


inlet pressure and temperature. The flow rate was 
initially calculated using: 

*7 Dfajrtfc -P Vv ) (12) 

However, a venturi for use in LOX service is 
typically calibrated in a room-temperature water 
facility 1 . This means that, due to the cryogenic 
temperature during motor firing, the actual venturi 
throat diameter is somewhat smaller - introducing a 
bias into the LOX flow rate measurements. One 
involved method for calculating the actual throat area 
would require the use of a finite element model. An 
easier method, however, is to simply use the 
coefficient of thermal expansion and the temperature 
difference between the calibration facility and the 
cryogenic feedline to reduce the throat diameter an 
appropriate amount, resulting in a slightly modified 
form of (12) 


( 10 ) 


The denominator of (10) results from 
accounting for all fluids that are heated and expelled 
through the motor nozzle. This includes, of course, 
the main LOX and fuel, but also includes the 
TEA/TEB ignition fluid and nitrogen purge masses. 

It should be noted that the time interval used 
to calculate each average factors out of equation (10). 
Additionally, the integrals in (10) are calculated using 
a trapezoidal rule approximation - except for the fuel 
total mass, which is calculated by direct fuel weights 
pre- and post-test. The appendix contains the general 
DRE used to calculate motor C*. 

Nearly all the terms in the DRE shown in the 
appendix can be broken down into constituent 
components. For example, A* is the simple average 
of the pre-and post-fire throat areas. But, to 
complicate things further, both of the diameters used 
in the average area calculation are themselves simple 
averages of four (each) pre- and post-fire 
measurements (see the equation below). 
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m*£[(l+tfJ')D r Yc d> jp(P V ' -P Ym ) 

(13) 

Equation (13) has two terms that are 
dependent upon the inlet temperature; namely, p and 
P Vlf . These two quantities are calculated using GASP 1 
and a Wagnerian curvefit 4 , respectively. ' 

Two types of ignition sources are used for 
the subscale motors. The 1 1-inch GOX motor uses a 
small GOX-hydrogen torch, and all three LOX 
motors use TEA/TEB, a pyrophoric fluid. Following 
a pre-set length of time over which the torch Is lit or 
TEA/TEB is flowing, a gaseous nitrogen purge is 
started to keep combustion products from flowing 
back into the ignition feed system. Therefore, (12) is 
used to calculate TEA/TEB flow rate, with the 
density modeled as a function of temperature only (a 
conceptual bias), and with the vapor pressure 
assumed to be zoo 5 . 

Nitrogen is metered using the same orifice as 
for the TEA/TEB (or the GOX-hydrogen torch). The 
equation for sonic gas flow through a sharp-edged 
orifice is 


The ramification of equation (11), explored 
in greater demil later in this paper, is that the throat 
area is assumed to change linearly - i.e., the average 
throat diameter changes as a second-order function of 
time. 

The LOX mass flow rate is metered using a 
cavi taring venturi, which requires the measurement of 
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The nitrogen y is also calculated using GASP. 

The fuel integral is a simple weight change 
calculation baaed upon average pre- and post-fire 
weighing: 




rK ] (15) 

'^+^+*0 


Differentiation of (15) shows clearly that the 
pre- and post-fire grain weight measurements have 
derivatives that are of opposite sign. This fact, 
highlighted in a previous section, helps to reduce the 
overall measurement uncertainty. The pre- and poet- 
fire measurements have correlated bias uncertainties, 
because the same load cell (or scale, for the lighter 
motor pieces) is used for both sets of weights. Since 
a bias uncertainty is systematic, the magnitude of the 
uncertainty does not change from pre- to post-test 
weightings. In effect, the bias is “calibrated out,” 
resulting in a more certain weight difference 
measurement than either individual pre- or post-fire 
weights. 

One additional point that should be noted is 
that the chamber pressure integral (and resulting 
average), was originally based upon measurements 
taken in the headend of the motor. A well-known 
phenomenon exists in solid and hybrid rocket motors, 
however, such that there can be a significant pressure 
drop down the length of the motor that is a complex 
function of chamber conditions and flow rate. With 
the addition of a mixing chamber aft of the fuel grain 
in hybrid motors (which promotes combustion 
product stagnation due to a diffuser- like effect and 
amplifies the pressure drop), it was felt that using 
headend pressure measurements contributed another 
conceptual bias to the performance calculations. 

Overall, the DRE for measured C* becomes 
a complicated amalgam of expressions (see appendix) 
that defies closed form partial derivative calculations. 

Theoretical C* 


investigated how measurement uncertainties affected 
the modeling of the SSME power balance. A similar 
process to that used in Coleman and Stem was used 
for the work reported in this paper, though no attempt 
was made to estimate the uncertainties in the 
theoretical calculations themselves. The uncertainty 
in the theoretical C*, then, results directly from 
uncertainties in the input data. 

Historically, theoretical C* was calculated 
by interpolating tabulated values that resulted Sum 
thermochemistry runs completed by a consortium 
member. The values were tabulated as a function of 
0/F and chamber pressure. However, the values did 
not account for the TEA/TEB, the GOX-hydrogen 
torch gases, the silica phenolic from the aft mixing 
chamber and the baffle (in the 11 -inch GOX motor) 
or the nitrogen purge. It was felt all these fluids 
contributed to the measured C* pressure integral and 
to the theoretical calculations, and that not accounting 
for them introduced additional biases into the 
experimental results. Therefore, this analysis used the 
latest release of the classic NASA-Lewis ODE 
module written by Gordon and McBride called CEA 1 , 
along with a complete input of all the propellants 
expended during a test 

Aft pressure measurements were, also used 
for input to CEA. Other inputs include motor 
propellant formulation — calculated during mixing 
and casting operations - and 0/F --. calculated using 
total expended propellants via direct weight 
measurements or integrated mass flow rates. 

Since a computer code is now used to 
calculate the theoretical C*, closed form solution for 
partial derivatives is rendered impossible. Therefore, 
they, too, were estimated numerically by perturbing 
Individual measurements which affected the CEA 
input data. 

C* Efficiency 

C* efficiency is a measure of how much 
energy is released by combustion - as gauged by the 
pressure generated - compared to the theoretical 
chemical potential energy stored in the pfopellantfs). 
A simple equation for combustion efficiency is 


Estimating the uncertainties associated with 
a theoretical model is a relatively new addition to the 
uncertainty analysis arena. Coleman and Stem 7 have 
investigated using uncertainty techniques to validate a 
CFD code in which the validation level is set by the 
combination of the uncertainties' in the experimental 
data and the portion of the uncertainties in the code 
that can be estimated Brown and Coleman 1 also 


Vr 


c* 

jj. ^ Measurtd 

~ c* 

^ ODE 


( 16 ) 


Since both the. denominator and the 
numerator are themselves functions of numerous 
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primitive variables and inputs, equation ( 16 ) is meant 
only to convey the notion of combustion efficiency. 

Computer Code 

I jsr Interfax 

Due to the requirement to integrate mass 
flow rates and the fact that closed form solutions for 
partial derivatives were rendered impossible, a 
computer program was written for this analysis. It is 
a Windows^-based, FORTRAN 90 routine that 
incorporates GASP and CEA. A Visual FORTRAN" 
GUI* was written to make the code as simple to use as 
possible. Figure 1 shows the opening screen for 
SHUE. 



Fig. 1: SHUE opening screen. 


SHUE will be capable of performing 28 
separate uncertainty analyses, ranging from mass 
flow rate and total impulse uncertainties to C*, I*, and 
Iv. (measured, theoretical and efficiency) 
uncertainties. Currently, all but the l n and ly m 
analyses are complete and validated. To date, there 
are approximately 25,000 SLOC, of which 8,000 are 
from GASP and CEA, Once completed, it is 
estimated SHUE will contain more than 30,000 
SLOC. 

The user is required to input which motor is 
being analyzed, which fUel formulation was used, and 
a test name, as well as pre- and post-test 
measurements for fuel grain weights and nozzle 
dimensions. Additional input includes fluid meter 
parameters, elemental bias limit estimates, and the 
data-sampling rate. A typical input dialog appears in 
Figure 2. 



Fig. 2: Typical SHUE dialog box. 


Once all of the required information is 
provided, the user selecis any or all of the analyses 
that are required to be performed This is selected 
through the drop-down menus at the top of the 
program window, an example of which is shown in 
Figure 3. 



Fig. 3 s Selection of an uncertainty analysis. 


After selecting which analyses to perform, 
the user then opens the data file for the motor test and 
selects “Gol” under the ‘Tile" menu. Depending 
upon the analysis, the user may be required to input 
■ start- and stop-times for the limits of. integration. 
Run-times vary from a few seconds to three hours. 
Analysis results are currently output to an ASCII text 
file. 

Algorithm Description 

The repetitive algorithm used for SHUE is 
straightforward, but because of the shear number of 
calculations, the algorithm is quite intricate. After the 
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user inputs the required data and times, the code loops 
through every measurement (pressures, temperatures, 
et cetera) to calculate a nominal value of the 
experimental result. Calculations for nozzle and 
weight measurements are completed during user 
input, and the integrals are calculated using the 
trapezoidal rule. Other, more sophisticated schemes 
for numerical integration were initially considered, 
but were abandoned because results for the 
trapezoidal rule were nearly identical and required 
less CPU time. 

Once the nominal result is calculated, the 
partial derivatives are estimated by singly perturbing 
every measurement used in the nominal calculation 
by 0.5% — hence the name "jitter” routine - and using 
equation (3). It should be noted that perturbing a 
single variable in an integral means that a complete 
integral is required with each pass through the loop. 
In other words, if there are 5000 pressure data points, 
a single point is perturbed and a new integral with 
that single perturbed variable (along with the other 
4999 nominal data points) is calculated. This feet 
required that the code be written using REAL*8 
variables in order to calculate the derivatives with 
sufficient accuracy. This jitter process was repeated 
for every variable that affected the experimental 
result. For example, there are 25 types of 
measurements which affect motor C* efficiency. Of 
those 25 types, 9 involve recalculating integrals for 
each of the multiple measurements taken in that type 
(venturi inlet pressure, for example), resulting in 
approximately 12,300 partial derivatives for a typical 
20-second test. 

Once each partial derivative is calculated the 
code passes control to the uncertainty propagation 
subroutine, where the partial derivatives are 
combined with the uncertainty estimates. Another 
parameter is also calculated in this portion of the 
code: the UPC. 

Most simply, the UPC allows the ranking of 
the most significant uncertainty contributors) by 
tabulating what percent of the overall uncertainty- 
squared each type of measurement contributes. 
Equation (17) shows how the UPC is calculated for 
each variable: 


UPC = 



*100 


(17) 


Results and Discussion 

24-Inch Motor 

In order initiate validation of $UUE, data 
from a 24-inch test (HP24-5030) was used for this 
analysis, os that test was a close analog to the 
upcoming 250K motor tests. Once SHUE is complete 
and validated, the HP24-5030 test data will be scaled 
to match expected conditions for the 250K test firing. 
The scaling will be done in order to get a pre-test 
uncertainty prediction, so that an understanding can 
be derived as to which measurements are most 
important for reducing the uncertainty in the full scale 
motor tests. 

There are 23 variables or measurements that 
are required to calculate C* (measured, theoretical, 
and efficiency). Factory specifications were used for 
estimating bias limits for pressure transducers, 
thermocouples, and load cells. Estimates for the 
venturi and orifice calibration biases were assumed 
based upon the results of an unpublished calibration 
analysis for the 250K motor venturi completed by 
engineers at MSEC. LOX p and nitrogen y biases 
were taken from the GASP user’s manual, as stated in 
[2]. The LOX vapor pressure bias estimate was taken 
as the largest reported error estimate in [4]. Finally, 
weight percent bias estimates for each constituent in 
the feel formulation was based upon a conservative 
estimate (Le., larger than is probable) of the bias in 
the scale used at the casting feefety divided by each . 
constituent’s total weight per mix. 

It should be noted that there are only three 
measurements that have the potential for significant 
uncertainties introduced from precision sources: the 
nozzle throat diameter measurements, the feel grain 
weights, and the headend weights. The nozzle throat 
measurements will be discussed separately, but 
typically, asymmetric erosion occurs during the motor 
firing. These asymmetries are treated like precision 
errors and combined with the bias limit for the 
calipers used, by using a root-sum-square method. 
Similar procedures are used for the grain and headend 
weights, though the precision contribution is small in 
both cases. Table 1 summarizes each measurement 
along with its nominal value and bias limit. 
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Parameter 

Nominal Value 

Bias Limit 


580 osla 

5 psl 

P 

1500 osia 

7 psi 

T 

-268 °P 

3°F 


1500 psia 

5 psi 

T 

71 T 

3*F 


65.5 PSi 

0,l%Readlne 

Pu* 

68 Ibm/ft* 

0 .15% Reading 


44 Ibm/ft 1 

10% Reading 

Tlttma 

1.4447 

0.5% Reading 

LOX Venturi C, 

0.9622 

0.019(2%) 

LOX Venturi a 

9.6*10'*/ < T 7 

10% Reading 

LOX Venturi 

7 

70 T 

5T 

N, Orifice C. 

0.30 

0.016(2%) 

Nj Orifice Diameter 

0.09 in. 

" 

0.005 in. 

Grain Weight Loss 

322 Ibtn 

25 Ibm' 

Headend Weight 
Loss 

23.4 lbra 

0,005 Ibm* 

D* 

4.6 In. 

0.010 in. 

Etcorez wtfV 

N/A 

lxl 0' 1 % 

R-45 wt% 

N/A 

2xlO J % 

DDI *t% 

N/A 

6x10'* % 

Voranol wt% 

N/A 

lxlQ-’W 

Agents White wt% 

N/A 

lxl0 J % 

Carbon Black. wt% 

N/A ‘ 

3x10’ % 


Table li Bias limit* for each experiment*! jjatameter. 


In general, the C* measurements are fairly 
good - C* efficiency had a ± 2 % uncertainty band. 
Table 2 summarizes the results of the analysis, and 
compares the old calculation method results with the 
new method proposed in this paper. The comparison 
shows that the changes in procedure affects the results 
at a magnitude lower than the overall uncertainty - 
i.e., the changes, for all practical purposes, cannot be 
seen for lack of measurement resolution* However, 
the “lessons learned” in this process have more 
significant ramifications for the 1 1-inch GOX motor. 



Old Vahie 

New Value 

Uncertainty 

c*_. 

5753 ft/s eo 

5775 ft/sec 

113 ft/sec 1 


5769 fVsac 

5772 fVsec 

14 ft/seo j 

*u 

99.7% 

100% 

2% 1 


Table 2; Overall uncertainty analysis results, 


The detailed breakdown of uncertainty 
contributions is interesting, and U probably more 
informative. Tables 3 through 5 summarize UPC 
values, with uncorrdated and correlated UPC 


1 These bias limits apply to each individual pro- and 
post-fire measurement and not to the overall weight 
difference (see Uncertainty Analysis for explanation). 
# Due to ITAR restrictions, these nominal values 
cannot be reported, 


components reported for all three experimental 
results. Note that six measurements (P^ w 

LOX Venturi ct, Pt*a> N 4 
Orifice C 4 ) do not appear in the following three 
tables. They each contribute less than 0.01% to the 
measurement uncertainty for the three experimental 
results. 


Parameter 

UPC 

Component^) 

Resultant 

Measurement 

UPC 

LOX Venturi C. 

51.04% " 

5L04% 

LOX Venturi T_ . 

2.12% 

2.12% 

Uncorrelated P 

0.01% ' 


Correlated 

19.07% 

19.08% 

Uncorrelated 

P— ' — 

0.00% 


Correlated 

0.86% 

0.86% 

p 



Uncorrtlated 

0.01% 


TflUCU-, 



Correlated 

5,44% 

• 5.45% 

T 



Unoorrelated p, w 

0.00% 


Correlated p lflK 

0.20% 

0.20% 

Uncompleted N, Orifice 

0,06% ’ 


Diameter 



Correlated N x Orifice 
Diameter 

0.07% 

0.13% 

Uncoirelated Pre-Fire 

0.20% 


D* 


Uncoirelated Post-Fire 

4.42% 


a* 



Correlated Pre-Fire D* 

. 0.60% 


Correlated Post-Fire D* 

13.27% 


Correlated Pre- and 

2.32% 

20,82% 

Post-Fire D* 



Uncorrelated Pre-Fire 

0.00% 


Headend Weight 



Uncorrelated Post-Fire 

0.00% 


Headend Weight 



Correlated Pre-Fire 

Headend Weight 

0,00% 


Correlated Poat-FIre 

0.00% 


Headend Weight 



Correlated Pre- and 

-0.00% 

0.00% 

Post-Fire Headend 

Weight 



Unccrrelated Pre-Fire 
Center Segment Weight 

36.71% 


Uncoirelated Post-Fire 

37.73% 


Center Segment Weight 



Correlated Pre-Fire 

73.42% 


Center Segment Weight 



Correlated Post-Fire 

75.47% 


Center Segment Weight 



Correlated Pre- and 

-223.10% 

0,24% 

Post-Fire Center 



Segment Weight 




Table 35 Detailed UPC value* for C* T — ^ 
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Correlated 


Uocontlated N, Orifice 
Diameter 


Correlated N, Orifice 
Diameter 


F 


Correlated Fuel 

Constituent 
Weight % 


Uncorrelated Pre-Fire 
Headend Weight 


Correlated Pre-Fire 
Headend Weight 


Correlated Post-Fire 
Headend Weight 


274.84% 


345.73% 


549.68% 


691.46% 


-1847,72% 


Table 4: Detailed UPC/UMF values tor 

It should be stressed that for all three cases, 
the most significant contributor to the overall 
experimental uncertainty is the LOX venturi C r In 
addition, the fuel grain and headend weight 
measurements are negligible contributors to the 
overall uncertainty. Some discussion provides some 
insight as to why those measurements once appeared 
to contribute significantly. 


CncorrcUtcd N, Orifice 
biameter 

0,10% 

Correlated N, Orifice 
biameter 

0.00% 

Uocontlated Fuel 

Constituent 
fright % 

0.02% 

Correlated . Fuel 

Constituent 

}Vei«ht% 

0.00% 

^correlated Pre-Fire 

0.24% 



Correlated Pro- and 
Post-Fire D* 



Correlated Poit-Fire 
Headend Weight 

0.00% 

Correlated Pro- and 
Poit-Fire Headend 

Weight 

-0.00% 


78.90% 


Uncorrelated Poit-Fire 86.43% 

Center Segment Weight ■ 

Correlated Pie-Fire 157.79% 

Center Segment Weight 

Correlated Poit-Fire 172.86% 

Center Segment Weight , . . 

Correlated Pro- and 494.98% 1.00% 

Poit-Fire Center 

Se pneru W eight - 

Table 5: Detailed UPC/UMF values for 

Inspection of the three tables shows a 
significant sensitivity to the overall result for either 
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pre- or post-fire weights - which seems intuitive. 
However, as explained in the Uncertainty Analysis 
section, the partial derivatives for the weight 
measurements are of opposite sign, resulting in the 
small UPC for all three experimental results, This is 
most dramatically evidenced by the enormous 
correlated bias contribution made by the. center 
segment weights. In effect, the weight difference can 
be measured to a much higher degree of accuracy 
than either individual pre- or post-fire weights. 

Another, significant bias contributor i$ the 
nozzle throat diameter measurements. The individual 
throat measurements are quite good: the caliper bias 
used for this analysis was conservatively estimated to 
be 0,010” - of which 0,001” comes from the caliper 
calibration, and the order of magnitude increase 
comes from estimated errors in being able to 
consistently measure the same point on the throat 
(depth, angle, etc.). The majority of that large overall 
contribution comes from the post- fire measurements. 

This is based upon the fact that nozzles with 
a significant erosion rate see a measurable 
asymmetric erosion pattern* Table 6 summarizes the 
pre- and post-fire measurements for the nozzle used 
in HP24-3030, As can be seen, there is significant 
variation in the measured diameter. 


Measurement 

Pre-Fire 

Post-Fire 

Clocking 

Diameter (in.) 

Diameter (in.) 

0*- 180 4 

3,706 

5.375 

45° -225° 

3.706 

5.278 

90° - 270 0 

3.706 

5.269 

135* -315 4 

3.706 

5.386 

Average 

. 3.706 

5.327 


Table <5: Nozzle throat measurements for HP24*5030. 


This diameter variation, in reality, represents 
a source of precision uncertainty. That is, until more 
quantitative data exist as to actual flowfield 
conditions the variance can be considered a random 
effect The : most common way to account for 
precision uncertainties is by standard statistical 
analysis, which was the method used for this analysis. 
While these aren’t strictly measurement uncertainties, 
they do fall into the conceptual bias category. This is 
because the C* equation assumes a circular throat 
area; in effect, forcing the calculated area to be 
circular despite the evidence to the contrary. Further, 
because the erosion patterns are asymmetric, the 
resulting average post-fire diameter will be 
completely dependent upon where the measurements 
are taken. For example, if the throat were 
inadvertently clocked by 30* and measured again, the 


average would be quite different. Therefore, the best 
way to estimate that bias uncertainty is by calculating 
the precision index, P x , (based upon the precision 
index of the mean), and treating that index as a 
fossilized bias. 

The overall post- fire throat diameter bias, 
then, was nearly an order of magnitude higher t hti fl 
the 0.010” bias from the actual measurement 
instrument* it was calculated as 0.099”. The authors 
assert that doubling or tripling the number of post-fire 
measurements would reduce the effects of this 
conceptual bias to a more acceptable level. 

11-Inch Motor Considerations 

Another issue that is related to nozzle throat 
erosion has stemmed from the results of some recent 
1 1-inch GOX motor firings. These motors appear to 
have impossibly high combustion efficiencies - 
107.4% for HP 11-4010, for example. As mentioned 
previously, using the aft pressure trace, accounting 
for the silica phenolic weight, and using CEA to 
calculate the theoretical C* reduces the efficiency to 
102.1%, This value is only 2% higher than the 
theoretical maximum, within the uncertainty band for 
the measured C+ calculation. However, the 11-inch 
GOX motors consistently run higher than 100% 
efficiencies, which suggests that there is a bias that 
has not been addressed Some have suggested that 
perhaps the average throat area cdculatioh needs 
adjusting for the QOX motors. In feet, the difficulty 
with the throat area calculation, historically, is the 
reason that some propulsion specialists prefer 1^ 
efficiency measurements 16 . 

The simplest change one could consider 
making to the average throat area calculation is to 
assume that throat diameter (not area) is a linear 
fraction of time. That is, assume the throat diameter 
increases at a constant, time- averaged rate. Doing so 
reduces the efficiency for HPl 1-4010 to 99,2%, a 
slightly more believable level of performance. 
However, it is difficult, considering the highly 
complex (and. mostly unexplored) heat transfer and 
erosion characteristics Qf hybrid motor nozzle throats, 
to assume a simple, linear phenomenon is controlling 
erosion. 

A different method to consider appears to be 
closer to the true nozzle response. A theory was 
postulated by one consortium member that assumes 
for a small, but measurable time, the nozzle wall acts 
much like a heat sink while it reaches the temperature 
required to ablate the graphite throat 11 - suggesting 
some parabolic or exponential function. To test this 


10 




TO WEEKS 


PAGE . 013/015 


'98 15:84 .-"CM LrCPCS PPCG OFC 


hypothesis, throat diameter histories for 3 recent 11- 
inch tests were curve-fit using an exponential function 
of the form 


d' 



(18) 


where a is equal to the initial throat diameter and b is 
a constant Parameters a and have direct physical 
significance, so to be a meaningful model, b should 
also have some physical significance. The actual 
throat conditions are difficult to measure, so four non- 
dimensional parameters were chosen to describe 
them, in anticipation that the parameter b would be a 
function of these parameters. 

First, to account for the fluid dynamic 
effects that contribute to mechanical erosive effects, 
the Reynolds number based upon diameter (Re^ was 
chosea Second, to account for the heat transfer 
effects, that contribute both to mechanical and 
chemical erosive effects, both the Nusselt number 
(where h, - based upon initial throat diameter - was 
estimated using the Bartz equation 1 *) and the Prandtl 
number were chosea And finally, because excess 
oxidizing species can enhance throat erosion due to 
the formation of CO and C0 2 through scavenging of 
the graphite, the sum of the mass fractions of O and 
Oj in the exhaust species was chosen. In order to 
factor Out the dependence of Re^ and Nu^ upon the 
throat diameter, the product of those two parameters 
was used (see the equation below). Additionally, all 
theoretical parameters (except the film coefficient, as 
described above) were taken from CEA output 


4 ” 


4A f m 


k jqjD m 


(19) 


Results of curvefitt between b and the non- 
dimensional parameters show a remarkable 
relationship. That final curvefit is shown in Equation 
( 20 ). 

b = <p Mn (NuKey- im Pr" 1 ” 11 (20) 

To generate confidence in the curvefit, test 
data from the JIRAD program were used to compare 
curvefit predictions from Eq, (20) with curvefits of 
throat diameter data using Eq, (18). Those results, 
summarized in Table 7, confirm that for this motor 
configuration, the constant b does have physical 
significance to actual motor conditions. This 


realization lends credibility to the claim that the throat 
erodes as an exponential function of time. It should 
be stressed that equation (20) is not yet capable of 
predicting throat erosion, as seen by the relatively 
high relative error for JIRAD 9. Additionally, 
another result for a test with an O/F near 
stoichiometric conditions showed virtually no 
agreement This effort was only to deduce whether 
(18) was an adequate engineering model, not a 
prediction tool. 


Motor 

b (Eq. 22) 

b (Eq. 24) 

% Error 

JIRAD 4 

47.7 

49.7 

4% 

JIRAD 9 

34.6 

46.7 

35% 

JIRAD 16 

32.3 

28.0 

13% 


Table 7: Comparison of erosion models. 


. As a final check, the Biot number (By the 
ratio between the thermal resistances to convection 
and conduction) for HP11-4G10 was calculated to be 
~700* which shows that the 1-D throat thermal 
gradient is very steep, farther indicating the throat 
indeed acts as a heat sink 

Table 8 shows how significantly the choice 
of models affects 


Linear Area 
Change 

Linear Diameter 
Change 

Exponential 
Diameter Change 

102.1% 

99.2% 

97.3%, 


Table 8: reusing the three erosion models. 


. These results suggest a need to begin using 
(18) for calculating average throat area. The erosion 
data have more than a coincidental relationship to die 
major non-dimensional parameters that govern heat 
transfer, Additionally, it has been suggested dial 
there is an unaccounted-for bias that most likely 
comes from throat area calculations. Therefore, the 
authors suggest using Equation (18) for all throat area 
calculations for the 1 1-inch motor. 

How this applies to the 24-inch and 250K 
motors is unclear, but will be investigated in the near 
future. A 3-D carbon cloth/phraplic throat is 
typically used for those motors, so its through-ply 
thermal conductivity will vary from that assumed for 
the graphite throat So, while the reliance on non- 
dimensional heat transfer and flow parameters will 
likely result, it is doubtful equation (20) is directly 
applicable. 
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Conclusions 

Measurement uncertainty was estimated for 
three performance parameters used in rocket 
propulsion testing. Results showed that measured C*, 
theoretical C*, and C* efficiency uncertainties are 
most sensitive to LOX mass flow rate, chamber 
pressure, and nozzle throat diameter. It was also 
shown that ftiel grain weight measurements are 
negligible contributors to these uncertainties. 
Additionally, increasing the number of post-fire 
nozzle throat measurements would reduce 
experimental uncertainties, which is a result of the 
significant asymmetric erosion typically found. 

Throat erosion models and throat area 
calculations were investigated as a possible source of 
conceptual biases which result in impossibly high C* 
efficiencies. The investigation showed that the nozzle 
throat in a hybrid motor, likely erodes as an 
exponential function of time, hot a linear function of 
time as previously assumed The erosion model was 
compared to a simple theoretical model consisting of 
four major heat transfer and fluid flow parameters, 
and showed good agreement These results indicate 
that that an exponential-type throat erosion model is 
fully based in physical parameters that govern heat 
transfer. 

. Future work on SHOE will continue, with 1^ 
and Iy K performance measurement capability added 
In addition, new throat erosion models will be 
incorporated, as well as improved bias source 
estimates and 250K motor configuration fidelity. 
Uncertainty estimates based upon expected operation 
for the 250K motor will also be generated Finally, 
throat erosion . modeling will continue, and will 
incorporate the more extensive database of 24-inch 
motor test data 
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